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^ , Abstract 

o : 

. Mirror matter is a possible dark matter candidate. It is predicted to exist if parity is 

D \ an unbroken symmetry of the vacuum. The existence of the mirror matter, which in 

Q^' addition to gravity is coupled to our world through photon-mirror photon mixing, 

1^ ■ would result in orthopositronium (o — Ps) to mirror orthopositronium (o — Ps') 

. . \ oscillations. The experimental signature of this effect is the invisible decay of o — Ps 

.5^ ^ m vacuum. 

^ \ This paper describes the design of the new experiment for a search for the o—Ps — > 

' invisible decay in vacuum with a sensitivity in the branching ratio of Br{o — Ps — > 

invisible) ~ 10"^, which is an order of magnitude better than the present limit 
on this decay mode from the Big Bang Nucleosynthesis. The experiment is based 
on a high-efficiency pulsed slow positron beam, which is also applicable for other 
experiments with o — Ps, and (with some modifications) for applied studies. Details 
of the experimental design and of a new pulsing method, as well as preliminary 
results on requirements for the pulsed beam components are presented. The effects 
of o — Ps collisions with the cavity walls as well as the influence of external fields 
on the o — Ps — > o — Ps' oscillation probability are also discussed. 

Key words: orthopositronium, invisible decay, mirror dark matter, pulse positron 
beam 
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1 Introduction 



The dark matter problem provides one of the strongest indications for physics 
beyond the Standard Model. Although the unknown physics is usually ad- 
dressed in a direct manner in high energy experiments, new results may also 
be expected from precision experiments at lower energies. 
Orthopositronium (o — Ps, the triplet e+e^-bound state), is a particularly 
interesting system for such an approach [1-3] . For example, it has been shown 
recently that experiments searching for invisible decays of o — Ps with the 
(currently achievable) level of sensitivity in the branching ratio Br{o — Ps 
invisible) ~ 10~^ — 10~^ have significant discovery potential [2]. An obser- 
vation of o — Ps —>■ invisible decay would unambiguously signal new physics 
phenomenon which could be induced either by the existence of extra dimen- 
sions [4], or of fractionally charged particles [5-7], or of hght gauge bosons 
[2]. Other interesting experiments with o — Ps are motivated by tests of high 
order QED corrections to the o — Ps decay rate [8], searching for a violation of 
fundamental symmetries in positronium annihilation [9], tests of antimatter 
gravity in the free gravitational fall of positronium [10], the possibility to ob- 
serve positronium Bose- Einstein condensation [11] and others. Among them 
the most exciting one is probably related to the search for the dark matter of 
a mirror-type. 

Mirror matter is predicted to exist if parity is an unbroken symmetry of nature. 
The idea was originally discussed by Lee and Yang [12] in 1956, who suggested 
that the transformation in the particle space corresponding to the space in- 
version X — > —X should not be the usual transformation P, but PR, where R 
corresponds the transformation of a particle (proton[12]) into a refiected state 
in the mirror particle space. After observation of parity nonconservation Lan- 
dau assumed [13] that R—C, i.e. he suggested to identify antiparticles with 
the mirror matter, but then CP must be conserved which we know is not 
the case. The idea was further developed by A. Salam [14], and was clearly 
formulated in 1966 as a concept of the Mirror Universe by Kobzarev, Okun 
and Pomeranchuk [15]. In this paper it was shown that ordinary and mirror 
matter can communicate predominantly through gravity and proposed that 
the mirror matter objects can be present in our universe. 

Since that time the mirror matter concept has found many interesting appli- 
cations and developments. In 1980 it has been boosted by superstring theories 
with X Eg symmetry, where the particles and the symmetry of interactions 
in each of the Eg groups are identical. Hence, the idea of mirror matter can 
be naturally accomodated. Today's mirror matter models exist in two basic 
versions. The symmetric version proposed early was further developed and put 
into a modern context by Foot, Lew and Volkas [16]. The asymmetric version 
was proposed by Berezhiani and Mohapatra [17]. More detailed discussions of 
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mirror matter models can be found in Ref.[18]. 



In the symmetric mirror model the idea is that for each ordinary particle, 
such as the photon, electron, proton and neutron, there is a corresponding 
mirror particle, of exactly the same mass as the ordinary particle. The PR 
operator interchanges the ordinary particles with the mirror particles so that 
the properties of the mirror particles completely mirror those of the ordinary 
particles. For example the mirror proton and mirror electron are stable and 
interact with the mirror photon in the same way in which the ordinary proton 
and electron interacts with the ordinary photons. The mirror particles are not 
produced in laboratory experiments just because they couple very weakly to 
the ordinary particles. In the modern language of gauge theories, the mirror 
particles are all singlets under the standard G = SU{3) ® SU{2)l ® U{1)y 
gauge interactions [16]. Instead the mirror particles interact with a set of 
mirror gauge particles, so that the gauge symmetry of the theory is doubled, 
i.e. G®G (the ordinary particles are, of course, singlets under the mirror gauge 
symmetry) [16]. Parity is conserved because the mirror particles experience 
^^ + ^4 (i.e. right-handed) mirror weak interactions while the ordinary particles 
experience the usual V — A (i.e. left-handed) weak interactions. 

It was realized some time ago by Glashow[19], that the orthopositronium sys- 
tem provides a sensitive way to search for the mirror matter. Glashow's idea 
is that if a small kinetic mixing of the ordinary and mirror photons exists 
[21], it would mix ordinary and mirror orthopositronium, leading to maximal 
orthopositronium - mirror orthopositronium oscillations, see Figure 1. Since 
mirror o — Ps' decays predominantly into three mirror photons these oscil- 
lations result in o — Ps invisible decay in vacuum. Remember, that due 
to the odd-parity under C transformation o — Ps decays predominantly into 
three photons. As compared to the singlet (1^5*0) state (parapositronium) , the 
small o — Ps decay rate (due to the phase-space and additional a suppression 
factors) gives an enhancement factor ~ 10^, making it more sensitive to an 
admixture of this new interaction [2,3,20]. 

Photon-mirror photon kinetic mixing is described by the interaction Lagrangian 
density 



where F^^ {F'^y) is the field strength tensor for electromagnetism (mirror elec- 
tromagnetism) . The effect of ordinary photon - mirror photon kinetic mixing 
is to give the mirror charged particles a small electric charge[16,19,21]. That 
is, they couple to ordinary photons with charge 2ee ^ . 

^ Note that the direct experimental bound on e from searches for 'milli-charged' 
particles is e ~ 10~^ [5,7]. 
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Fig. 1. The double degeneracy between orthopositronium mass eigenstates of ordi- 
nary o-Ps and mirror o-Ps ' is broken when a small mixing ( upper picture ) term is 
included. 

Orthopositronium is connected via a one-photon annihilation diagram to its 
mirror version (o — Ps'). This breaks the degeneracy between o — Ps and 
o — Ps' so that the vacuum energy eigenstates are (o — Ps + o — Ps') / \/2 
and (o — Ps — o — Ps')/\/2, which are spht in energy by AE = 2hef, where 
/ = 8.7 X 10^ MHz is the contribution to the ortho-para sphtting from the 
one-photon annihilation diagram involving o — Ps [19] . Thus, the interaction 
eigenstates are maximal combinations of mass eigenstates which implies that 
o — Ps oscillates into o — Ps' with probability: 

P(o - Ps ^ o - Ps') = sin^ out, (2) 

where cu = 27re/. 

In the simplest case of o — Ps — > o — Ps' oscillations in vacuum[19] , because 
the mirror decays are not detected, this leads to an apparent increase in the 
decay rate, since the number of o — Ps, satisfies 

N = cos^ a;te~^^'"* ~ exp[-t{rsm + uj^t)], (3) 

where Tsm = 7.039934(10)/is^"'^ [22] is the Standard Model decay rate of o — 
Ps (i.e. when the oscillation length goes to infinity). Thus T^^^-^ ~ ^smij- + 

The above calculation is not applicable to an experiment, where the positron- 

ium is confined in a cavity, because in this case the collision rate is not zero and 
the loss of coherence due to the collisions must be included in the calculation 
of Eq.(3) [23,24]. Note, that the probability P{o - Ps ^ o - Ps') can also be 
affected by an additional splitting of o — Ps and o — Ps' states by an external 
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electric or magnetic field [23]. This is quite similar to the phenomenon of n — n 
[25] or muonium to antimuonium oscillations [26] in various environments, see 
section 2 for a discussion of these effects. 

Recently, Foot discussed implications of the DAMA [27] and CRESST [28] 
experiments for mirror matter-type dark matter which is coupled to ordinary 
matter through the interaction of Eq.(l) [29] (for references related to mir- 
ror matter see e.g. [30]). He has shown that the annual modulation signal 
obtained by the DAMA/Nal experiment, as well as the CRESST data, can 
be explained by mirror matter-type dark matter if the photon-mirror photon 
mixing strength is in the region 

e ~ 4 X 10-^ (4) 



Interestingly, this value of e is also consistent with all other known experimen- 
tal and cosmological bounds including SN1987a ^ and the standard Big Bang 
Nucleosynthesis (BBN) bound [33]. It is also in the range of naturally small 
e- values motivated by grand unification models [18]. 

If e is as large as in Eq.(4), the branching ratio Br{o — Ps invisible) for 
invisible decays of orthopositronium in vacuum is of the order (see Section 4): 

Br{o - Ps ^ invisible) ~ 2 x 10~^ (5) 

For comparison, the BBN limits [33] deduced from the successful prediction 
of the primordial ^He abundance are 

e < 3 X 10-^ (6) 

and 

Br{o - Ps ^ invisible) < 10"^ (7) 

respectively. 

The first experiment on the o — Ps — > invisible decay, motivated by a puzzle 
in the o — Ps decay rate (see below), was performed a long time ago [34], and 
then repeated with higher sensitivity [35]. The results exclude contributions 
to the o — Ps decay rate from invisible decay modes (such as o — Ps — >• uu, 
miUicharged particles, etc..) at the level of Br{o — Ps — > invisible) < 3 • 10~^, 
but are not very sensitive to the o — Ps — > o — Ps' oscillation mechanism 
because of the high collision rate in these experiments. Indeed the limit on 
e extracted from the results of rcf. [35], taking into account the suppression 
collision factor, is e < 10~^ [23] and is not strong enough to exclude a possi- 
ble mirror-matter contribution at the level comparable with the BBN limit of 

^ The SN1987a limit e < 10~^'^ obtained in Ref.[31] is actually much weaker. For 
a more detailed discussion of this constraint, see Ref.[32] 
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Eq.(6). 

Given the indications for the mirror world coming from dark matter [29] and 
the neutrino physics [36,37], as well as the intuitive expectation that nature 
could be left-right symmetric, it is obviously important to determine exper- 
imentally whether orthopositronium is a window on the mirror world. Since 
there arc no firm predictions for e. experimental searches for o— Ps invisible 
have to be performed with a sensitivity as high as possible. In particular, the 
question whether the sensitivity for the interesting branching ratio of Eq.(5) 
is experimentally reachable has to be studied. 

The goal of this work is to present a design of the experiment for a search 
for o — Ps invisible decay in vacuum with a sensitivity better than the 
corresponding BBN limit Br{o — Ps invisible) < 10~^ and high enough 
to check the prediction of Ref.[29], Br{o — Ps invisible) ~ 2 x 10~^. The 
experiment requires the production and subsequent decays of o — Ps's to oc- 
cur in vacuum [23,24], and hence, the use of a specially designed slow positron 
beam operating preferably in a pulsed mode to enhance the signal-to-noise 
ratio for the efficient tagging of o — Ps production [38]. The experimental 
signature of o — Ps — > invisible decay is the absence of an energy deposition 
which is expected from the ordinary o-Ps annihilation in a 47r calorimeter 
surrounding the positronium formation target. We show that this signature is 
clean and that the signal events can be identified with a high confidence level 
due to the efficient tagging of the positron appearance in the target and the 
high-efficiency measurement of its annihilation energy. It should be noted that 
recent substantial efforts devoted to the theoretical and experimental deter- 
mination of the QED o — Ps properties, hopefully results in a solution of the 
long-standing discrepancy between the measured and predicted orthopositro- 
nium decay rate in vacuum, see e.g. [39,40,41]. However, the current level of 
the theoretical precision achieved by Adkins ct al. [8,22] is about two orders 
of magnitude better than the experimental one [40,41,42,43]. Thus, further 
positron beam based experiments to measure the o — Ps decay rate in vacuum 
are required and are of great interest to test high-order QED corrections. 
The paper is organized as follows. In section 2 we consider the o—Ps o—Ps' 
oscillations and effects of various environments on it. In section 3 we report on 
the design of the experimental setup to search for the o—Ps — > invisible decay 
in vacuum. The description of the detector components including the design 
of a high-efficiency pulsed positron beam is presented in sections 3 and 4. The 
preliminary simulation results and expected sensitivity level are discussed in 
sections 5 and 6, respectively, section 7 contains concluding remarks. 
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2 o — Ps ^ invisible decay rate in a vacuum cavity 



In the simplest case of o — Ps — > o — Ps' oscillations in vacuum [19] the branch- 
ing ratio occurring during a long enough observation time can be calculated 
as 



Br{o -Ps^ invisible) = (S) 

Eq. (8) may not be applicable to an experiment in a cavity. It is well known that 
collisions damp the oscillations, e.g. in the limit where the collision rate is much 
larger than the decay rate (or oscillation frequency, whichever is smaller) the 
effect of the oscillations becomes negligible. In addition, external fields might 
result in a loss of coherence due to additional splitting of mass eigenstates. 
Thus, their effect must be included [23,24]. 

Let us first consider the case where the collision rate is much larger than the 
decay rate, Tcou ^ [24], then the evolution of the number of orthopositro- 
nium states, N, satisfies: 

^ ~ -TsmN - T^uNp, (9) 

where the second term is the rate at which o-Ps oscillates into o — Ps' (whose 
subsequent decays are not detected). In this term, p denotes the average os- 
cillation probability over the collision time. That is. 



P 



rcoii f'e-^-^^''sm'u;t'dt' ^ r,„, fe-"^--'' {ivt'fdf , (10) 
Jo Jo 



where we have used the constraint that the oscillation probability is small, i.e. 
cut <^ 1. As long as t ^ l/F coii, a reasonable approximation for the vacuum 
cavity experiment, then 

P-^. (11) 

coll 

Thus, substituting the above equation into Eq.(9) we have 

F'ff ~ r,^ + = r,^ 1 + . (12) 

coll \ coll'- sm ) 

The difference between the higher decay rate measured in the vacuum cavity 
experiment, relative to the value predicted by theory, can be expressed as 

r -r — ri3^ 

^ exp ^ sm — -p p V / 

coll'- sm 
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For the cavity size and the o — Ps emission spectrum in a recent experiment on 
the o — Ps decay rate in vacuum [42], we estimate that Vcoii ^ ^^oPs-, which, 
neglecting the contribution from external fields (which are in fact negligible 
in this case[23]), implies that 

u;^ ~ 2 X lO-^r^p^ e < 10"^ (14) 

Thus, the hmit of Eq.(14) is still not strong enough compared to the BBN one 
of Eq.(6). 

In the presence of a static external electromagnetic field the o — Ps — o — Ps' 
degeneracy is broken and the probability that positronium decays as mirror 
o — Ps', rather than ordinary o — Ps, is 

Br{o -Ps^ invisible) ^ ^4^27//)^ + ^''^ 

where A represents an additional oscillations damping factor combining effects 
of collisions of o — Ps's with the cavity walls, scattering on residual gas atoms 
and the influence of external flelds [44] . 

Estimates obtained taking into account Zeemann and Stark effects in positro- 
nium and o — Ps scattering in the Van der Waals potential of residual gas 
moleculas in the cavity show that for the external magnetic held ~ 100 G, 
the electric held ~ 100 V and residual vacuum pressure < 10~^ Torr, result 
in a value A < Fsm [44]. The most crucial parameter is the number N of 
o — Ps collisions during its lifetime in the cavity, which gives a suppresion 
factor ~ 1/iV. Hence should be as small as possible, i.e. N < 1. 



3 Experimental setup to search for the o — Ps ^ invisible decay 

The experiment is designed with the goal to observe the o — Ps invisible 
decays, if its branching ratio is greater than 10~^. Figure 2 shows a schematic 
view of the experimental setup in which positrons from a pulsed beam [38] are 
stopped in the MgO target and either form positronium, i.e. o — Ps or p — Ps, 
or annihilate promptly into 27's. The secondary electrons (SE) produced by 
the positrons hitting the target are accelerated by the voltage applied to the 
target relative the grounded transport tube. Then they are transported by a 
magnetic held in the backward direction relative to the positrons moving in 
spirals along the magnetic field lines and defiected to a microchannel plate 
(MCP) by a, E X B filter as is shown in Figure 3. 

The trigger for data acquisition is generated by a coincidence within ±3 ns 
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Fig. 2. Schematic diagram of the experimental setup. 

of a pulse from the MCP and the signal from the pulsed beam, which is 
synchronized with the positron arrival time at the target. 

Accordingly, the apparatus is designed with several distinct and separated 
parts: i) a pulsed slow positron beam and a low-mass target for efficient or- 
thopositronium production in a vacuum cavity, ii) a positron appearance tag- 
ging system with a high signal-to-noise ratio based on a high performance 
MCP, iii) an almost Arc BGO crystal calorimeter (ECAL) surrounding the 
vacuum cavity for efficient detection of annihilation photons. The cavity has 
as little wall mass as possible to minimize photon energy absorption. 

The occurrence of the o — Ps —>■ o — Ps' — > invisible conversion would appear 
as an excess of events with energy deposition comparable with zero in the 
calorimeter above those expected from the prediction of the background. In 
case of a signal observation the number of excess events could be cross-checked 
by small variations of experimental conditions which affect the o—Ps — > o—Ps' 
transition rate but do not result in a loss of energy from ordinary positron an- 
nihilations. The identification of signal events relies on a high-efficiency mea- 
surement of the energy deposition from the annihilation of positrons. Com- 
pared to the previous experiments on o — Ps — > invisible decay , the use of the 
BGO calorimeter gives a significant reduction in total required ECAL mass 
(see section 3.3). 
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To achieve a sensitivity in the branching ratio of 10~^ in a reasonable amount 
of data-taking time, the rate of o — Ps decays per second has to be as high as 
possible consistent with minimal reduction of the o — Ps —>■ invisible signal 
efficiency and acceptably small dead time. For the pulsed positron beam design 
presented in Section 4, the trigger rate in the photon detector is expected to 
be ~ 100 Hz which is low enough to allow these events to be recorded without 
losses. 

3.1 Positron tagging system 

The SEs produced by positrons hitting the a — Ps production target are used 
to tag the time of positron appearance in the target. The positron tagging 
system is based on a high performance MCP as a SE detector. The low energy 
SE emerging from the target are accelerated by an electric field and deflected 
to the MCP by the E x B fllter, as shown in Figure 2. The system works in 
detail as follows. The pulsed beam of positrons with energy ~ 500 eV is guided 
by a magnetic field with a value of S ~ 100 G and passes through a region 
with crossed electric and magnetic fields {E x B region 1 in the schematic 
diagram of Figure 2). The transversal electric field value is ~ 500 V/cm. 
Positrons drift in the crossed electric and magnetic fields with a velocity given 
by 

Vd^Ex B/B'' (16) 

For the given values of the electric and magnetic fields the drift velocity is 
y^^TxlO^m/s resulting in the positron displacement about of 11 mm in 
the drift region 1. The electric field in region 2 has the same value but the 
opposite direction relative to region 1. As a result, downstream of the region 
2 positrons will move back to the axis if the transport system, separated from 
slow electrons and ions. Then, positrons are transported to the target in the 
curved magnetic field created by the coils. Preliminary results illustrating the 
method are shown in Figure 3, where the calculated trajectories of secondary 
electrons passing trough the E x B filter are shown in the Y-Z plane for a 
wide primary positron beam. The energy spectrum of secondaries is taken 
from Ref.[45], the angular distribution is assumed to be isotropic. 

The target is a disk with a diameter of ~ 10 mm and a thickness of 0.1 mm. 
The SE acceleration potential of -300 V is applied to the target. The tar- 
get surface is coated by MgO for an efficient production of orthopositronium 
and an efficient secondary electron emission [45] . The secondary electrons are 
transported in the backward direction relative to the positrons, see Figure 
3, moving in spirals along the magnetic field lines. Thus, trajectories of the 
secondary electrons will be spatially separated from the positron trajectories 
by a distance estimated to be ~ 25 mm, significantly larger then the diameter 
of the positron and the test electron beams (~ 5 mm). The MCP installed 
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between region 1 and region 2 (see Figure 2) detects the electrons. The back- 
ground count rate of the MCP can hmit the efficiency of the positron tagging. 
This and other sources of the background can be suppressed by an appropriate 
choice of the MCP type and by using a pulsed positron beam with a low duty 
cycle. The energy spectrum and the yield of SE are important for the tagging 




Fig. 3. Calculated trajectories of secondary electrons in the tagging system in Y-Z 
plane. The primary positron beam is increased in size for illustration purpose. 

efficiency. The SE emission is a surface effect, involving only a very thin layer 
of the target material in the process. Thus, the SE yield and the MCP out- 
put signal are proportional to the energy loss dE/dx of positron. The SE are 
emitted with an energy up to 100 eV and a large angular spread. This results 
in a dispersion of their arrival time at the MCP. To reduce this effect the SE 
should be accelerated immediately after production by the electric field at the 
biased target. Note that the MCP pulse shape could be used for additional 
discrimination between the signals from positrons and the noise pulses. 

The factors affecting the performance of the positron tagging system are 

• the coefficient of secondary electron emission; 

• the efficiency of the electron transport from the target to the MCP, and the 
efficiency of the MCP itself; 

• the MCP noise level and the environmental background; 

• physical backgrounds, from e.g. beam interactions with residual gas and 
cavity walls, with material of the Ex B ffiters etc. accompanied by electron 
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or ion production. 

We have selected the MCP Hamamatsu F4655-12 for the tagging system be- 
cause it provides the best signal-to- noise ratio. The MCP signal rise time is 
~ 300 ps. The signal amplitude varies typically from 0.1 to 1 V depending 

on the applied voltage, the MCP gain and the number of the emitted SEs. 
The collection voltage between the MCP bottom surface and the anode is not 
critical, a setting above 100 V is sufficient for > 95% collection efficiency. 

3.2 o — Ps production target 

A MgO coated film is planned to be used as o — Ps production target. The 
production rate of o — Ps per single positron is about 25 % at the positron 
energy of ~ 700 eV [41]. The MgO is also reported to have the high secondary 
electron emission coefficient ( > 4 ) at this energy. Another possibility is a 
high porosity and low density Si02 target [42]. 

To ensure a minimal probability of total absorption of annihilation photons, 
the target itself and the surrounding components of the target region must be 
carefully designed. To minimize the thickness of dead material, in particular 
the thickness of the beam pipe in this region is very important. Based on 
simulation results a 1 mm thick Al pipe seems the best choice ( see section 
5). One possible way to avoid the problem of dead material is to drill a hole 
in the central crystal (Figure 2) and glue the beam pipe into it in such a way 
that the a — Ps production and decay region is surrounded almost completely 
by the crystal active medium . 

3.3 Photon detector 

The optimal choice of the 7-dctcctor (ECAL) can be made by the following 
considerations. The total ECAL mass W is given roughly by 

W ~ An/SpL^ (17) 

where p and L are respectively mass density and the radius of the ECAL 
detector. We chose L ~ 2OA511, where A511 is the attenuation length of 511 keV 
7's. The relevant parameters for different types of materials used in ECAL's 
are listed below in Table 1. 

The required mass is minimal for a BGO ECAL due to its high effective Z 
(remember, the photo- absorption cross-section a ~ Z^). Another important 
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Table 1 

Comparison between different types of ECAL. 



ECAL type 


BGO 


Nal 


CsI(Tl) 


Scintillator 
plastic/liquid 


A511 


~1 cm 


~2.5 cm 


~1.9 cm 


~10 cm 


density, g/crrfi 


7.1 


3.6 


4.5 


1.0 


ECAL mass, kg 


~240 


-189 


~1034 


-33510 


A^/511 keV 


~ 4 • 10^ 


~ 20 • 10^ 


~ 10 • 10^ 


~ 2 • 10^ 


Hygroscopic 


no 


yes 


slightly 


no 



feature of BGO's is that they are not hygroscopic, thus, no additional dead 
material has to be introduced. 



The schematic drawing of the 7-detector is shown in Figure 2. The 7-quanta 
produced in positron or positronium annihilation are detected by a (almost) 
47r BGO crystal calorimeter [46] . We plan to use BGO crystals also for mea- 
suring of the photon time with respect to the arrival time, of Ih^ positron 
bunch on the target. The detector system consists of about 100 BGO crystal 
surrounding the vacuum beam pipe as shown in Figure 2. The full system is 
calibrated and monitored internally using the 511 keV annihilation line. For 
the crystal wrapped in aluminized mylar the light yield was measured to be 
200±14 photoclectrons/1 McV. This results in a probability of zero energy 
detection due to Poisson fluctuation of the number of photoelectrons, to be 
less than 10~^^ for the zero energy signal defined as events with energy depo- 
sition less than 50 keV [47,48] . This result justifies the selection of the BGO as 
the 7-detector. The crystals, which have been lent to us by the Paul Scherrer 
Institute (Villigcn, Switzerland), have a hexagonal shape with a length of 20 
cm and an outer diameter of 5.5 cm, their original wrapping is a 0.75 mm 
thick tefion. In order to reduce this amount of dead material, the inner ring 
of BGO's has been wrapped in a 2 //m thick foil, aluminized from both side 
with 1000 A thick layers. The required number of crystals, determined with 
the simulations, provides an almost isotropically uniform thickness of 20-22 
cm of BGO. 



4 High-eflficiency pulsed positron beam 

In this section we report the preliminary design of a high-efficiency pulsed 
slow positron beam for particle physics experiments with orthopositronium 
in vacuum. Our primary consideration is that the system should be of the 
magnetic transport type because this provides the simplest way to transport a 
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slow positron beam from the positron source to its target [51]. An electrostatic 
beam is also considered, but presently the problems with construction, time 
schedule and increased cost dominate over the benefits. The basic working 
principals involved in the design of magnetic or electrostatic variable-energy 
DC positron beams are well known [52]. The advantages and disadvantages 
of a beam formed and transported by a magnetic field in comparison to the 
electrostatic one are also known, see e.g. [53]. 

Various techniques to create pulsed positron beams have been reported with 
the main focus so far on material science applications [50]. Those attained 
by the Munich [54] and the Tsukuba [55] groups use RF power in the pulsed 
beam formation. However, acceleration by the sine function of the RF electric 
fields is by no means the optimal choice. The system based on this method 
requires a wide time window of chopping and accordingly the beam efficiency 
becomes small. 

A pulsing system with a higher performance has recently been proposed by 
Oshima et al.[56]. The main idea is the same as for the RF method: to adjust 
the time-of-flight for each positron according to the time it arrives at the 
accelerating point. However, instead of applying a sinusoidal RF field, a more 
suitable pulse shape of the electric field is generated such as an approximate 
inverse parabol function of time [56] . This method has been further developed 
by lijima et al [57] for the material measurements in which the lifetime or 
time-of-flight of orthopositronium atoms is close to the o — Ps lifetime in 
vacuum ~ 142 ns. For these applications it is necessary to modify the originally 
proposed technique [56] in order to generate higher intensity positron beams 
by accumulating positrons over a wider time span even though the bunch 
width becomes larger, but still much less than the typical measured timing 
intervals of ~ 100 ns. By using a high permeability buncher core a bunch width 
of 2.2 ns (FWHM) for 50 ns collection time and a repetition period of 960 ns 
has been achieved [57] . The main problem encountered in this technique is the 
limitation of the voltage supplied by a post-amplifier to the buncher. 

Our new pulsing method allows to compress an initial positron pulse of 300 ns 
to about 2 ns pulse width. The method relies on a positron velocity modulation 
by using a new double gap buncher technique. The detailed description of the 
system components as well as their requirements are presented. 



4.1 Design criteria and overall system design 

The beam is supposed to be used for several different experiments with o — 
Ps in vacuum. Thus, the final beam construction should compromise several 
design goals which are crutial for them and which are summarized as follows: 
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• simple and not very costly experimental apparatus, 

• beam energy range from 100 eV to 1000 eV, 

• beam intensity of ~ 10^ — 10^ positrons per second, 

• pulse duration at the target Stx < 3 ns for an initial pulse duration at the 
moderator 6tM — 300 — 400 ns, 

• repetition rate 0.3-1.0 MHz, 

• high peak/noise ratio, (single) Gaussian shape of the pulse, 

• beam spot size at the target position is of the order of a few millimeters 
assuming 3-5 mm ^^Na source diameter, 

• minimal pumping time of the vacuum system. 




Fig. 4. Schematic illustration of the magnetically transported pulsed positron beam. 

Figure 4 shows schematic illustration of the pulsed positron beam design. The 
positron pulsing section consists of a chopper and a buncher and is based 
on positron velocity modulation combined with the RF bunching technique. 
A positive potential (~ 100 eV) is applied to the moderator foil in order 
to insure the proper energy of the positrons at the buncher input. Initial 
positron pulses with duration 300 ns are formed with the chopper grid placed 
2 mm apart from the moderator foil. The pulsed voltage with an amplitude 
about of +5V applied to the chopper grid relative the moderator foil will 
stop slow positrons with energy about 3 eV emitted from the moderator. Fast 
positrons emitted from the source are eliminated from the beam by the velocity 
analyzer (90 degrees curved solenoid, placed downstream the chopper). When 
the voltage applied to the chopper grid is zero, the positrons come through 
the chopper grid and are accelerated in the gap between the chopper grid and 
first drift tube (see Figure 4). Thus, positron pulses with a duration of 300 ns 
are produced by this way. The chopper voltage pulses can be produced by a 
standard fast-signal generator. 
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Fig. 5. Schematic illustration of the W(IOO) single crystal moderator. 

In the gap between the drift tube 1 and the buncher tube the velocity of 
positrons from the 300 ns pulse is modulated by a nonlinear pulsed voltage 
applied to the buncher tube relative to the drift tubes. The buncher tube 
length is determined by a distance-of-flight of positrons entering the buncher 
during 300 ns. In a second gap between the buncher tube and a drift tube 2 
the positron velocity is modulated again by the same voltage pulse applied to 
the buncher. 

The buncher voltage pulse is produced by an arbitrary waveform generator 
(AWG) connected to a fast post-amplifier. The pulse shape for the two-gap 
buncher is determined by calculations described in Section 3. In accordance 
with Liouville's theorem the compression ratio in this case is determined by 
the ratio of the final and initial energy dispersion in the positron beam pulse. 
Experimentally measured initial energy dispersion of the moderated positrons 
is about 2 eV (FW at 10% of maximum) [61]. Taking into account that the 
final energy spread in the given two gap buncher is about 200 eV we get the 
expected compression ratio of ~ 100. 

The moderation and associated voltages applied as shown in Figure 5 pro- 
duce a positron beam which contains both moderated positrons with energy 
of about 3 eV and an unmoderated {j3^ decay energy) positron component 
which has to be rejected. There are several possibilities to build a positron 
velocity filter to select the slow, moderated component and reject high energy 
component from the source. We consider at the moment a simple bending 
filter, which is made of solenoidal coils positioned one after the other. Com- 
pensation of the beam drift due to the curved magnetic field is also foreseen. 

The emittance of the beam is defined by the area in the phase space (r, a) 
multiplied by the square root of the energy and divided by a factor vr, where r 
and a are the radial positron position and angle of the beam. The brightness 
of the beam B can be defined as [52] B = ,2^2^ ; where / is the intensity, E 
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the particle energy, the angular divergence of the beam, and D the beam 
diameter. The brightness of the beam is limited by its initial emittance and by 
Liouville's theorem. This means that if the beam size is decreased, the angular 
divergence is increased. One of the problems is that positrons derived from 
the moderator will gain in emittance because of the presence of a transverse 
electric field in the vicinity of the moderator mesh, see Figure 5. 



4.2 Analytical design of the pulsing system 

To avoid the already mentioned over-voltage problem with the post-amplifier, 
we try first to design a simple pulsing system which accelerates and/or decel- 
erates positrons only in a single gap during 300 ns, such that positrons arrive 
at the target after the buncher almost simultaneously. This problem can be 
solved exactly if all positrons at the entrance of the buncher have exactly the 
same longitudinal momentum. We use a time dependent electric field in the 
gap so that particles arriving first are decelerated and the later ones are accel- 
erated. For a given distance from the gap to the target we have in this case one 
free parameter - the deceleration potential at the time t = 0. There is a soft 
limitation in the choice of the initial value of the deceleration potential: the 
momentum of positrons coming first should not be very small, since a small 
smearing of the initial momentum would cause a large dispersion in time at 
the target. However, we found the accelerating gap potential near the end of 
the 300 ns time interval should be more than 1 kV, which is difficult to achieve 
in a high frequency pulsed mode operation. 

More convenient and economic is the buncher with two modulation gaps. In 
this case, as is shown in Figure 4, the buncher consists of the entrance drift tube 
1, buncher electrode and exit drift tube 2. The decelerating or accelerating 
potential is applied to the buncher electrode with respect to the drift tubes 
resulting in positron velocity modulation in two gaps. 

The simulations of the extraction optics, beam transportation and of the ve- 
locity modulation of positrons are performed with the GEANT4 [63] and 3D- 
Bfield [62] codes with the goal to minimize the timing resolution and optimize 
the shape of the bunching pulse. 

The numerical solution for the shape of the buncher pulse is shown in Figure 
6. It was calculated for the following characteristics of the system: 

• moderated positrons have an average kinetic energy of ~ 3 ±0.5 eV and are 
emitted isotropically from the fiat film of 10 mm in diameter. According to 
the previous measurements [61] the low energy tail with intensity of 10% 
has been added. 
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Fig. 6. The bunching voltages seen by positrons at the first and the second velocity 
modulation gaps, respectively. 

• the velocity filter tube axis is curved with a radius of 45 cm. Particles enter 
into the tube along Y (upward) and exit along the Z axis, 

• the buncher gaps are both 1 mm wide with the electric field along the Z 
axis. 

• the buncher electrode is 140 cm long, 

• the beam transportation tube is 1.05 m long and 10 cm in diameter 

The optimal duration of the buncher pulse and the shape of the potential 
which positrons "see" in the gaps are chosen to satisfy the following criteria: 

• the time difference between two positrons arriving at the gaps should be 
smaller at the second gap for the applied potential, 

• the amplitude of the RF pulse should be within ±60 V, 

• after modulation at the second gap, positrons should arrive at the target at 
the same time. 

We chose the parabolic time-dependent potential proportional to and chang- 
ing from from -60 V (decelerating part) to 30 V for the positron velocity mod- 
ulation at the first gap. Once this is fixed, the only free parameter left is the 
distance from the second gap to the target. The time dependence of the po- 
tential at the central electrode of the buncher at time t > 300 ns can then be 
calculated to make particles arriving at the target simultaneously. It is also 
assumed that the potential at the electrode at the end of the bunching pulse 
returns to its initial value -60 V. The pulse shape at the second gap is found 
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Fig. 7. Distribution of the positron time- of- flight at the target for the following cuts 
on the longitudinal energy E\ | of positrons emitted from the moderator: a) > eV, 
b) >1.5 eV c) >2.0 eV and d)> 2.5 eV. The initial positron energy distribution is 
Gaussian with an average of 3 eV and a = 0.5 eV, the angular distribution is taken 
to be isotropic. 

as the result of an iteration procedure for the solution of the corresponding 

equations. Figure 6 shows the resulting shape of the bunching voltages seen by 
positrons at the first and the second velocity modulation gaps, respectively. 



4.3 Requirements for the system components 



The time shape of the positron pulse at the target is affected by the perfor- 
mance of the pulsed beam components. In this section requirements for the 
monochromaticity of the moderated positrons and for the amplifier used to 
generate the buncher pulse shown in Figure 6 are considered. 



4.3.1 Monochromaticity of the moderated positrons 

In Figure 7 the simulated positron time distribution at the target is shown for 
different energy spread of the initial positrons on the moderator. The degree 
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Fig. 8. Residual response of the post amplifier 100A250A for the input pulse shown 
in Figure 3. 

of monochromaticity of the initial positrons was defined by the cuts on the 
longitudinal positron kinetic energy E\\. The angular distribution of initial 
positrons was taken for simplicity to be isotropic. The results show that the 
best resolution achieved with this method is about 2 ns. The results also 
show that the quality of the moderator is an important parameter. So, the 
annealing of the W-foil in situ is important. Otherwise quick degradation of 
the surface quality through interactions with a gas results in a practically 
isotropic reemission of the positrons, and hence a significant increase in the 
phase space and timing spread. 



4-3.2 Buncher pulse and amplifier responce 

The post-amplifier choice is made on the basis of the following requirements. 
We need an amplification factor of about 100 since the typical amplitude of 
an AWG (arbitrary waveform generator) is 1 V and the output signal has a 
peak to peak value of 90 V. The gain has to be programmable for final tuning 
of the output amplitude. The amplifier has to be sufficiently powerful to drive 
the 50 Ohm load. A few tens of Watts is expected for a typical signal rate. A 
wide frequency band from a few kHz to a few 100 MHz is required to minimize 
the signal shape distortions. And finally, the integral non-linearity has to be 
within at most 1% . 

We select commercial RF post-amplifier 100A250A (Amplifier Research [64]) 
as a supplier of a fast pulse to the buncher. It has 100 W output power and 
a frequency bandwidth from 10 kHz to 250 MHz. It is designed for 50 Ohm 
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loading with reliable over-current and over- voltage protection. 

The 100A250A post-amplifier is an AC-coupled device. To shift the signal 
basehne to the -60 V level, an external DC source is used. This voltage supplier 
is also used for fine pedestal tuning in order to compensate the baseline shift 
in the case of a high signal rate. 

The response of 100A250A was simulated following the circuit characteristics. 
The output pulse shape differs from the input due to the finite frequency 

bandwidth. The residual shape is defined as R{t) = Sout{t) — Sin{t) where 
Sin{t) is the input signal supplied by AWG with unit amplitude and Soutit) 
is the amplifier output pulse calculated for unit gain. Figure 8 shows R{t) 
expressed in percents of the input computed for the signal shape shown in 
Figure 6. The input shape is parameterized by an analytical function which 
has singularities (infinite first derivative) at t = 300 ns and t = 500 ns. Due 
to the upper frequency limit, the output signal is smoothed and differs from 
the input shape at these time points. This is indicated by two peaks on the 
residual plot. The post-signal baseline shift of about 1% is due to the lower 
frequency limit. 

One can see that the overall deviation of the response is not more than about 
±1%. It is expected that this value will be better for real signals because 
the AWG can not reproduce the singularities of the theoretical shape. The 
simulations of the beam showed that 1% of signal deviation does not result 
in a significant distortion of the bunched positron pulse shape. The RMS of 
the corresponding distribution (see Figure 7, d) has been changed by less than 
2%. However, for a deviation of the order of 5%, the RMS degrades from 2.3 
to 2.8 ns. This means that the shape of the buncher pulse must be reproduced 
within about ±1% of the theoretical shape. This value seemes to be achievable 
[64]. 



5 Monte Carlo simulations 

The positron trajectories in the beam were simulated with the GEANT4 and 
the 3D-Bfield programs [62]. The o— Ps production, propagation in the beam 
pipe, refiection on pipe walls and decay was simulated. The events for the 
o — Ps — > 37 process were generated taking into account the decay matrix 
element and assuming decays at rest. 

The Monte Carlo simulation of the photon detection in the apparatus was 
based on the GEANT3 package [63]. The geometries of the beam transport 
pipe, photon detector, positron tagging system and its material were coded 
into simulations. The simulation results were also benchmarked with the re- 
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suits of our previous experiment on a search for the o-Ps— > 7 + Xi + X2 decay 
mode [48]. In Figure 9, the distributions of the total energy deposited in 0.3 
mm thick stainless steel pipe (upper plot) and in a 1 mm thick Al pipe are 
shown for the 27 annihilation events. The 7-detection inefficiency was found to 
be less than 3 x 10~^ for the case of an Al pipe assuming that the zero-energy 
signal is defined as an event with the total energy deposition in the ECAL, 
E < 100 keV. 
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Fig. 9. Distributions of the energy deposited in a 0.3 mm thick stainless steel pipe 
(upper plot) and in a 1 mm, thick Al pipe. The total number of simulated, 2j-events 
is 10^ in both cases. The peaks at 511 keV and 1022 keV correspond to the total 
photo- absorption either of a single 511 keV photon or of both of them, respectively. 

Simulations show that the main contribution to the 7-inefficiency comes from 
the total (due to photo-absorption) or fractional (due to Compton effect) 
photon energy loss in the material of the vacuum beam pipe. In the latter 
case, the energy deposition in the ECAL corresponds to a small amount of 
photoelectrons, so that the ECAL signal could be comparable with "zero- 
energy" due to statistical fluctuations. The 7-detection inefficiency due to the 
limited ECAL thickness is estimated to be less than 10~^. A special mechanical 
design discussed in section 3.2 is necessary to reduce absorption in the pipe 
by choosing the appropriate materials and construction of the target region. 

The energy spectrum of o — Ps's produced at the MgO target is taken from 
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ref.[49]. The angular distribution is assumed to be isotropic and the intensity 
of the fast backscattered component [42] is taken to be 1% with a flat energy 
spectrum from 10 to 100 eV for a primary positron energy of a few keV. The 
coUisional o — Ps dissociation probabihty for a o — Ps kinetic energy greater 
than 10 eV is taken to be 100%. The average number of o — Ps coUisions 
with the cavity walls during their lifetime in the cavity was estimated to be 
Ncou ~ 2.5, resulting in corresponding suppression factor in the sensitivity 
of the experiment (see sections 2,6). The other concerns are about o — Ps 
decays near the entrance aperture to the detector, i.e. in the region of lower 7- 
detection efficiency, and the disappearence effect, namely when o — Ps's escape 
the cavity region through the entrance aperture without being detected. 



6 Sensitivity 



The experimental signature of the o — Ps invisible decay is an excess of 
events above the background at zero-energy deposition in the ECAL. The 
90%-confidence level limit on the branching ratio for the o — Ps — > invisible 
deacy for a background free experiment is given by 

^, ^, Nio — Ps ^ invisible) , , 

5(90%) = ^ i (18) 

^^o-Ps^^coll 

where N{o — Ps invisible) = 2.3 and the terms in the denominator are the 
integrated number of produced o — Ps's (No-Ps), and the average number of 
o — Ps collisions in the cavity, respectively. The number A^o-Ps is defined as 
a product Aq-Ps — Re+ " ^o-Ps " • t-i where the first factor is the number of 
delivered positrons per second on the target, the second one is the efficiency for 
o — Ps production, and the third one is the efficiency of the secondary electron 
transportation from the target to the MCP in the positron tagging system. 
Taking =2x W^/sec, Eo-ps = 20% and ee+ = 100%, we expect f« 7 x 10^ 
prompt and 1.7 x 10^ o — Ps annihilations per day. Thus, 5(90%) ~ 10~^. 

Eq.(18) gives the sensitivity for a background free experiment. We expect 
backgrounds which originates from the following sources: i) a fake positron 
tagging ii) the annihilation energy loss and iii) the disappearence effect. The 
are several sources of background that simulate the positron appearance signal: 

• the MCP after-pulses or pulses produced due to ionization of the residual 
gas by passing positrons; 

• cosmic rays; 

• environmental radioactivity. 

The MCP noise is typically low. Its level depends, e.g. on the "^^K contamina- 
tion in the surrounding materials and on the intensity of penetrating photons. 
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although the sensitivity of the MCP to photons is quite low. However, for 
the purpose of this experiment the S/N ratio must be better than ~ 10^. To 
reach this level a high SE emission coefficient of the target combined with a 
good MCP energy resolution and a high o — Ps production rate is crucial. 
A preliminary estimate shows that taking into account the noise spectrum of 
MCP4655-12 and the average number of secondary electrons per positron for 
the MgO target to be Nse — 4 [45], the MCP S/N ratio can be expected to 
be better than 10^. An additional noise suppression factor of ~100 (at least) 
is expected from the use of the MCP signal in coincidence with the positron 
beam pulse [38] with the time resolution of ~ 3 ns. The main source of beam 
associated background is expected from electrons and ions due to ionization 
of the residual gas atoms by positrons. Thus, a good vacuum is important. 
We estimate that with i?e+ — positrons on target, ~ 4 secondaries per 

positron and a vacuum in the cavity of ~ 10~^ Torr, a MCP S/N ratio > 10^ 
is achievable. The contribution from disappearence effect was also found to 
be small, ^ 5 x 10~*, although not completely free of assumption on fraction 
of fast o — Ps^s formed at the target. The preliminary overall background 
estimate results in Rbckg — 1-2 events /day. 

In case of the observation of zero-energy events, one of the approaches would 
be to measure their number as a function of the residual gas pressure in the 
cavity. This would allow a good cross-check: relatively small variations of gas 
pressure results in larger peak variations at zero energy due to the damping 
of o — Ps — > o — Ps' oscillations. 



7 Summciry 

The design of our experiment to search for the o — Ps ^ invisible decay 
mode in vacuum has been presented. The sensitivity of the experiment to the 
branching ratio of the o — Ps — > invisible decay in vacuum around ~ 10~^ 
seems to be achievable. Assuming that o — Ps ^ o — Ps' oscillations occur 
with mixing strength value e ~ 4 x 10~^ a total number of ~ 100 signal events 
would be expected in the ECAL during one month of data taking. Given that 
about 40 background events are expected a discovery with significance ~ 11 
could be possible [65]. 

One of the main features of the experiment is the use of a high-efficiency pulsed 
positron beam. The new proposed pulsing method allows to compress a 300 
ns initial positron pulse into a pulse with < 3 ns width. This will allow to 
enhance the signal-to-noise ratio for the efficient tagging of o — Ps production 
by more than an order of magnitude. A modification of the pulsing system 
could be made to use the beam for applications for material studies. Different 
components of the detector have to be constructed and fully simulated, in 
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order to test several crucial points of the design. 
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